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Abstract. Unit two-variable-per-inequality (UTVPI) constraints form 
one of the largest class of integer constraints which are polynomial time 
solvable (unless P=NP). There is considerable interest in their use for 
constraint solving, abstract interpretation, spatial databases, and theo- 
rem proving. In this paper we develop a new incremental algorithm for 
UTVPI constraint satisfaction and implication checking that requires 
0{m + n log n+ p) time and 0{n + m + p) space to incrementally check 
satisfiability of m UTVPI constraints on n variables and check implica- 
tion of p UTVPI constraints. 



1 Introduction 

The unit two-variable-per-inequality (UTVPI) constraints form one of the largest 
class of integer constraints which are polynomial time solvable (unless P=NP). 
There is considerable interest in their use for constraint solving [7,6], abstract 
interpretation [9], spatial databases [11] and theorem proving [8]. In this paper 
we develop new incremental algorithms for UTVPI constraint satisfaction and 
implication. 

A UTVPI constraint has the form ax+by < d where x, y are integer variables, 
c? S Z and a, 6 G {—1, 0, 1}. For example x + y < 2, x~y < —1, < — 1 and x < 2 
are UTVPI constraints. UTVPI constraint solving is based on transitive closure: 
A constraint ax — y < di and y + bz < ^2 implies the constraint ax + bz < di+d2- 
We can determine all the UTVPI consequences of a set of UTVPI constraints by 
transitive closure, but we need to tighten some constraints. The transitive closure 
procedure can generate constraints of the form x + x < d and —x ~ x < d, which 
need to be tightened to a; < [|J and —x < [|J respectively. 

JafFar et al. [7] and Harvey et al. [6] present incremental consistency checking 
algorithms for adding a UTVPI constraint c to a set (j) of UTVPI constraints. 
They are based on maintaining the transitive and tight closure of the set of 
UTVPI constraints 4> involving n variables. Both algorithms require ©(n^) time 
and 0{n^) space for an incremental satisfaction check. Both algorithms can also 
be used to incrementally check implication of UTVPI constraints by <^ U {c}. 
These algorithms require 0{n^ +p) time and 0{n'^ +p) space for an incremental 
implication checking, where p is the number of constraints we need to check 



for implication. In order to (non-incrementally) check satisfiability of m UTVPI 
constraints on n variables these approaches require 0{'n?m) time, and to check 
implication they require 0{n'^m + p) time. 

An improvement on the complexity of (non-incremental) satisfiability for 
UTVPI constraints was devised by Lahiri and Musuvathi [8]. They define a 
non-incremental satisfiability algorithm requiring 0{nm) time and 0{n -\- m) 
space. The key behind their approach is to map UTVPI constraints to difference 
constraints (also called separation theory constraints) of the form x — y < d, 
where x and y are integer variables and d G Z. 

The difference constraints are a well studied class of constraints because 
of their connection to shortest path problems. We can consider the constraint 
X — y < d as a directed edge x —>■ y with weight d. Satisfiability of difference 
constraints corresponds to the problem of negative weight cycle detection, and 
implication of difference constraints corresponds to finding shortest paths (see 
e.g. [3] for details). 

The mapping of UTVPI to difference constraints by Lahiri and Musuvathi [8] 
is a relaxation of the problem. The relaxed problem is solved by a negative 
(weight) cycle detection algorithm but it guarantees only the satisfiability in 
Q for the UTVPI problem. In order to check satisfiability in Z they need to 
construct an auxiliary graph and check for certain paths in this graph. 

In this paper we first extend Lahiri and Musuvathi's algorithm [8] to check 
satisfaction incrementally in 0{nlogn + m) and 0{n + m) space. Then we show 
how to build an incremental satisfiability and implication algorithm using the 
relaxation of Lahiri and Musuvathi and incremental approaches to implication 
for difference constraints of Cotton and Maler [3] , which can incrementally check 
implication in 0{n log n + m+p) time and 0{n + m + p) space. 

2 Preliminaries 

In this section we given notation and preliminary concepts. 

A weighted directed graph G = {V, E) is made up vertices V and a set E of 

weighted directed edges (u, u, d) from vertex u & V to vertex v & V with weight 

d. We also use the notation u —> v to denote the edge {u, v, d). 

A path P from vq to in graph G, denoted Vq v^, is a sequence of 
edges ei, . . . , efc where = {vi-i,Vi, di) G E. A simple path P is a path where 
Vi Vj,0 < i < j < k. 

A (simple) cycle P is a path P where Vq = and Wj ^ Vj,0 < i < j A fc A (i ^ 
OVjVfc). 

The path weight of a path P, denoted w(p) is E^^^di. 

Let G be a graph without negative weight cycles, that is without a cycle P 
where w{P) < 0. Then we can define the shortest path from vo to Vk, which we 
denote by SP{vo,Vk), as the (simple) path P from vo to Vk such that w{P) is 
minimized. 

Let wSP{x, y) = w{SP{x, y)) or +oo if no path exists from x to y. 



Given a graph G and vertex x define the functions ,5^ : F ^ M as 

6^ (y) = wSP{y, x) and (y) ^ wSP{x,y) . 

Let G be a graph without negative weight cycles. Then tt is a valid potential 
function for G if n{u) + d — ^{v) > for every edge (w, v, d) in G. 

There are many algorithms (sec e.g. [2]) for detecting negative weight cycles 
in a weighted directed graph, which either detect a cycle or determine a valid 
potential function for the graph. 

Given a valid potential function tt for graph G = {V, E) we can define the 
reduced cost graph rc{G) as (V, {(x, y, Tr{x)+d—Tr{y) \ {x, y, d) G E}). All weights 
in the reduced cost graph are non-negative and we can recover the original path 
length w{P) for path P from x to y from paths in the reduced cost graph since 
w{P) = w + 7r(y) — 7r(a;) where w is the weight of the corresponding path in the 
reduced cost graph. 

Since edges in the reduced cost graph arc non- negative wc can use Dijkstra's 
algorithm to calculate the shortest paths in the reduced cost graph in time 
0(nlogn + m) instead of 0{nm). 



2.1 Difference constraints 

Difference constraints have the form x—y<d where x and y are integer variables 
and d G Z. We can map difference constraints to a weighted directed graph. 

Definition 1. Let C be a set of difference constraints and let G = (V, E) he the 

graph comprised of one weighted edge x y for every constraint x — y < d in 
C. We call G the constraint graph of C. 

The following well-known result characterizes how the constraint graph can 
be used for satisfiability and implication checking of difference constraints. 

Tiieorem 1 ([3]). Let C be a set of difference constraints and G its corre- 
sponding graph. C is satisfiable iff G has no negative weight cycles, and if C is 
satisfiable then C\=x — y<diff w5P(a;, y) < d. 



2.2 UTVPI constraints 



A UTVPI constraint is of the form ax + by < d, where x and y are integer 
variables, a, 6 e {—1, 0, 1}, c e {—1, 1} and G Z. 

Definition 2. The transitive closure TC{4>) of a set of UTVPI constraints <p is 
defined as the smallest set S containing (j) such that 



ax 



cy < di G S Acy + bz < d2 G S 



ax + bz < di + d2 G S 



The tightened closure TI{(j>) of a set of UTVPI constraints <p is defined as the 
smallest set S containing (j) such that 



ax + ax < d G S 



ax < 



gS, a €{-1,1} 



Table 1. Transformation from UTVPI constraint c to associated difference constraints 
D(c) to edges in the constraint graph E{c). 



UTVPI c 


Diff. Constr. D(c) 


Edges E{c) 


X — y < d 


x'^ — y^ < d 
y~ — x~ < d 


id j_ 

y X 

_ d _ 

X ^y 


X + y < d 


x'^ — y~ < d 
y'^ — x^ < d 


y x^ 

- d + 

X ^y^ 


—X — y < d 


x~ — y'^ < d 
y—x^<d 


+ d _ 
y^ ^ X 

+ d _ 
a; -» 1/ 


X < d 


x+ -X- < 2d 


_ 2d + 

X ^ x^ 


—X < d 


x~ -x+ < 2d 


+ 2d - 
X^ ^ X 



The tightened transitive closure TTC{^) of (j) is the smallest set containing cp 
that satisfies both conditions. 

The fundamental result for UTVPI constraints solving is: 

Theorem 2 ([7]). Let <p be a set of UTVPI constraints. Then (j) is unsatisfiable 
iff exists 0<de TTC{<j)) where < 0. 

We can extend this for implication checking straightforwardly: 

Corollary 1. Let (j) be a satisfiable set of UTVPI constraints. Then (j) \= ax + 
by <d iff ax + by <d' & TTC{(i}) with d' < d or {ax < di,by < d2} C TTC{(i}) 
with dx+ d2 < d. 

Example 1. Consider the UTVPI constraints = {x — y < 2, x + y < — 1, 
—X — z < —4}, Then TC{(j)) includes in addition {x + x < 1, —y — z < —2, 
y — z < —5, ~z — z < —l^x — z < —3}. And TI{TC{(f>)) includes in addition 
{x < 0,-z < -4} and TTC{4)) = TI{TC((j))) in this case. The constraint 
—z < —3 is implied by </> as is y — 2: < 0. 

3 Lahiri and Musuvathi's approach 

Lahiri and Musuvathi map UTVPI constraints </> to difference constraints or 
equivalently a weighted directed graph G^, and they use graph algorithms to 
detect satisfiability. 

We denote the constraint graph arising from (p as = {V, E). The graph G 

contains two vertices x'^ and x~ for every variable x. These variables arc used to 
convert UTVPI constraints into difference constraints. The vertex x'^ represents 
+x and x~ represents —x. 

Let (/) be a set of UTVPI constraints. Each UTVPI constraint c ^ 4>\s mapped 
to a set of difference constraints -D(c), or equivalently a set of weighted edges 
E{c). The mapping is shown in the Table 1. Each UTVPI constraint on two 
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Fig. 1. (a) G^i for cj)' of Example 2 which is Q feasible but not Z feasible, (b) a zero 
length cycle in G^i. (c) for <p of Example 1. 

variables generates two difference constraints and accordingly two edges in the 
constraint graph. Each UTVPI constraint on a single variable generates a single 
constraint, and hence a single edge. 

Let —V denote the counterpart of a vertex v € V, i.e. — a;+ := x~ and 

—x~ := a;+. Clearly, for each edge {x,y,d) G E the graph G<{, also includes the 
edge (— 2/, —x, d) with equal weight. This correspondence extends to paths. 

Lemma 1 ([8]). // there is a path P from u to v in the constraint graph G^, 
then there is a path P' from, —v to —u such that w{P) = w{P'). 

If we relax the restriction on variables to take values in Z and allow them to 
take values in Q we can check satisfiability in Q using G^. 

Lemma 2 ([8]). ^ set of UTVPI constraints 4> is unsatisfiable in Q if and only 
if the constraint graph G^ = {V, E) contains a negative weight cycle. 

The reason why the feasibility in Z cannot be tested with arises from the 
possible implication of constraints of the form x + x < d or —x — x < d through 
the transitivity of constraints in <j). If d is odd (equivalently d/2 € Q \ Z) then (j) 
may be satisfiable with x = d/2 but not with x = [d/2\. 

Example 2. Consider the UTVPI problem 0' = {x — y < 2, x + y < —1, —x — z < 
—4, —X + z < 3}, then a transitive consequence of the first two is a; + a; < 1, 
while a consequence of the second two is —x — x < —1. Together these require 
X = 2- 

The graph G,j,' is shown in Figure 1(a). A zero length cycle is extracted in 
Figure 1(b). This cycle has solutions in Q but not in Z. □ 

The satisfiability algorithm of Lahiri and Musuvathi [8] is based on Lemma 2 
and the following result. 

Lemma 3 ([8]). Suppose G^f, has no negative cycles and (j) is unsatisfiable in 
Z. Then contains a zero weight cycle containing vertices u and —u such that 
wSP{u, ~u) is odd. 

The algorithm first checks Q feasibility using a negative cycle detection al- 
gorithm, and then checks that no such zero weight cycles exists in G^. 




Algorithm 1: LaMu 




Input: (f) a set of UTVPI constraints 




Output: SAT if 4) is satisfiable, UNSAT otherwise 


1 


Construct the constraint graph G4, — {V, E) from cj)\ 


2 


Run a negative cycle detection algorithm on G^; 


3 


if contains a negative cycle then 


4 


return UNSAT 


5 


else 


6 


let TT be a valid potential function for G,^ 


7 


E' := {{u, v) 1 (w, V, d) € E, 'n{u) + d = 7r(t;)}; 


8 


G; := {V,E')- 


g 


Group the vertices in G'^ into strongly connected components (SCCs). Vertices 




u and V are in the same SCC if and only if there is a path from u to v and a 




path from 1; to u in G'^. u and v are in the same SCC exactly when there is a 




zero- weight cycle in G4, containing u and v.; 


10 


for aW u £V do 


11 


if —u is in the same SCC as u and 7r(— w) — Tt{u) is odd then 


12 


return UNSAT 


13 


return SAT 



Example 3. A valid potential function for the graph shown in Figure 1(a) is 
7r(2/+) = 0, 7r(a;+) = 2, 7r(z+) = 5, 77(2/") = 3, 7r(a;-) = 1, 7r(z-) = -2. Each of 
the arcs is tight, so E' contains ah edges, and aU nodes are in the same SCC. 
Both a;+ and x~ occur in the same SCC and SP{x'^ ,x~) = 'k{x~) — 'k{x'^) = — 1 
is odd, hence the system is unsatisfiable. 

The complexity is 0{nm) time and 0(n + m) space assuming we use a 
Bellman- Ford single source shortest path algorithm [1,4] for negative cycle de- 
tection. 

4 Incremental UTVPI Satisfaction 

The incremental satisfiability problem is: Given a satisfiable set of UTVPI (j) 
(with 1} variables and m, constraints) and UTVPI constraint c, determine if 0U{c} 
is satisfiable. In this section we define an incremental satisfiability checker for 
UTVPI constraints that requires 0{n\ogn+m) time and 0{n+m) space. It relies 
on simply making incremental the algorithm LaMu of Lahiri and Musuvathi. 

The key is to incrementalize the negative cycle detection. We use an algorithm 
due to Frigioni et al. [5], using the simplified form (Algorithm 2: IncConDiff) 
of Cotton and Maler [3] (since we are not interested in edge deletion). Given 

a graph G = {V,E) and valid potential function tt for G and edge e = u -i 
V, this algorithm returns G' = (V, S U {e}) and a valid potential function tt' 
for G' or determines a negative cycle and returns UNSAT. The complexity is 
0(nlogn + m) time and 0(n + m) space using Fibonacci heaps to implement 
argmin. 



Algorithm 2: IncConDiff 

Input: = (V, E) a graph, tt a valid potential function for G<^, edge (m, d) a 

new constraint to add to G^. 
Output: UNSAT ii <p\J {u — V < d} is unsatisfiable, or G0u{u-t)<<i} and a valid 
potential function tt' for G^u{u-u<d}- 

1 7(v) := 7r(w) + d = 7r(t;); 

2 7(w) := for all w ^ v; 

3 while min{'y) < A 7(w) = do 

4 s := argmin(7) ; 

5 7r'(s) := 7r(s) +7(s) ; 

6 j{s) ~ ; 

7 for all s ^ t € G do 

8 if 7r'{t) = nit) then 

9 jit) := min{7(t), 7r'(s) + d' - n'{t)} 

10 if 7('u) < then 

11 return UNSAT 

12 return ((V.EU {(«,?;, d)}),7r') 



The incremental UTVPI satisfiability algorithm simply runs IncConDiff 
at step 2 of LaMu, the remainder of the algorithm is unchanged. Since the 
remainder of the the LaMu algorithm requires 0{n + m) time and space, the 
complexity bounds are the same as for the incremental negative cycle detection 
algorithm. 

5 Incremental UTVPI Implication 

The incremental implication problem is given a set P of p UTVPI constraints and 
a satisfiable set (j)oim UTVPI constraints on n variables, where <j) ^ c',yc' e P, 
as well as a single new UTVPI constraint c, check for each c' G P if (p Ac\= c'. 

Incremental implication is important if wc wish to use UTPVI constraints in 
a Satisfiability Modulo Theories (SMT) solver [10], as well as for uses in abstract 
interpretation and spatial databases. Our approach to incremental implication 
is similar to the approach of Cotton and Maler [3] for incremental implication 
for difference constraints. 

The key to the algorithm are the following three results. 

Lemma 4. Let ax + by < d £ TTC{(f)) where {a, b] C {—1, 1}, then ax + by < 
d&TC{(t>). 

The result holds since tightening introduces constraints involving a single 
variable and any further transitive closure involving them can only create new 
constraints involving a single variable. 

Lemma 5. Let ax<de TTC{(I)) where a e {-1, 1} then ax < d G TI{TC{4>)). 



The result holds since any result of transitive closure on a new UTVPI con- 
straint by < d! introduced by tightening, can be mimicked using the constraint 
by + by < {2d', 2d' + 1} that introduced it, and tightening the end result. 

The above two results show that TC{(j)) is the crucial set of interest for 
UTVPI implication checking. The following result shows how we can use the 
constraint graph to reason about TC{<j)). 

Lemma 6. c G TC{(j)) iff there is a cycle of length d, in the case of c = < d, 
or a path u v of length d in where {u, v, d) £ E{c). 

Proof. This lemma follows straightforward of the definition of TC{d>) and the 
transformation of the set of UTVPI constraints to its constraint graph. It can be 
prove easily by induction over the number of transitive closure steps in TC{(j)) 
resp. the length of cycle or path in G^. □ 

Example 4- Consider (j) of Example 1. Then for example x + x < 1 G TC{(j)) 
and there is a path x~ a;"*" of length 1 in shown in Figure 1(c). Similarly 
y — z < —5 G TC{^) and there are paths y~ ~^ z~ and ~^ of length —5 
in G^. 

We can use paths (in particular shortest paths) in to reason about most 
constraints in TTC{4'). In order to handle tightening we introduce a hounds 
function p which records the upper and lower bounds for each variable x, on the 
vertices a;+ and x~. It is defined as: 

wSP{u, —u) 
2 

We can show that p{x~) computes the the upper bound of x and — p(a;+) is 
the lower bound of x. Using Lemmas 5 and 6 we have. 

Lemma 7. For UTVPI constraints (j), 

p{x~) = min{d \ x<de TTC{(j))} 
p{x+) = min{d I -x<dG TTG{<t>)} 

where we assume min0 = +00. 

Example 5. Consider the graph in Figure 1(c) for constraints (p of Example 1. 
Then p{x~) = since wSP{x~,x+) equals to 1 and x < € TTG{^), while 

p{z+) = -4 since wSP{z+,z-) = -7 and -z < -4 £ TTC{(t>). Note e.g. 
p(x^) — +00 and there is no constraint of the form —x < d in TTC{(j)). 

The key to incremental satisfaction is the following result. 

Theorem 3. If the constraint graph G^ contains no negative weight cycle (i.e. 
<p is satisfiable in Q) then (p is unsatisfiable in Z iff a vertex v ^ V exists with 
p{v)+p{-v) < 0. 



p{u) = 



Proof. Let (j) he a, satisfiable set of UTVPI constraints in Q. Because of the 
Lemma 6 it applies the non-existence of a constraint < d € TC{4>) where d < 0. 
Therefore (j) is unsatisfiable in Z iff such a constraint belongs to TTC{(j))\TC{(j)) 
(Theorem 2), i.e. a possible unsatisfiability is caused by tightening. 

The Lemma 4 implies the equivalence for each constraint c G TTC{(j))\TC{4>) 
to ax < d where a G {—1, 0, 1}. Hence, (f> is unsatisfiable in Z iff two constraints 
X < di and —x < d2 with di + d2 < exist in TTC{(j)) iff (Lemma 7) p{x^) + 
p{x-) < 0. □ 

Effectively failure can only be caused by tightening if the bounds of a single 
variable contradict. 

Example 6. Consider the graph in Figure 1(a) for constraints 0' of Example 2. 
There is no negative weight cycle in G^i but p{x^) = and p(x+) = —1 because 

of x"*" — > ^ . Hence the system is unsatisfiable. 

Similarly the key to incremental implication is the following rephrasing of 
Corollary 1. 

Theorem 4. If ip is a satisfiable set of UTVPI constraints then 4> \= c iff for all 
{u, V, d) e E{c) either wSP{u, v) <d or p{u) + p{—v) < d. 

Proof. Let be a satisfiable set of UTVPI constraints. Because of Corollary 1 
it holds (j) ^ c and c = ax + by < d iS ax + by < d' G TTC{4)) and d' < d ov 
{ax <di,by< da} C TTC{<t)) and di + da < d. 

Now, the theorem holds straightforward due to Lemma 7 for the constraints 
with one variable, and Lemma 4 and 6 for the other constraints. □ 

Example 7. Consider the graph in Figure 1(c) for constraints (j) of Example 1. 
1= — 2; < — 3 is shown since iuSP{z~^, z~) = — 7 < 2 x —3. 

Algorithm 3 shows the new algorithm. As input it takes the constraint graph 
G^, a valid potential function tt, the bounds function p, a set P of UTVPI 
constraints to check for implication, as well as the UTVPI constraint c which 
should be added to (p. 

In the first step (line 1) the constraint c is transformed to its corresponding 
edges E{c) in a constraint graph. Then each edge in E{c) is added consecutively 
to the constraint graph C?^ by using the IncConDiff algorithm of Cotton and 
Maler [3]. After inserting all edges in G', the constraint graph equals to G^ujc} 
and tt' is its valid potential function for G' . Hence i^U {c} is satisfiable in Q. The 
remainder of the algorithm maintains the bounds function p' (lines from 8 to 12) 
and it is used to test the feasibility in Z (lines 13 and 14), and the implication 
of constraints in P (lines 15 to 18). 

By Lemma 6 to maintain p we need to sec if the shortest path from x to —x 
has changed. We only need to scan for new shortest paths using the newly added 
edges. We can restrict attention to a single added edge (u, v, d) since if there is 

a path from x over the edge (m, v, d) to —x (a;+ ^ u v x~) then because of 



Algorithm 3: ScSt - Incremental satisfiability and implication for UTVPI 
constraints. 




Input: = (V, E) a constraint graph representing set of UTVPI constraints 0, 




TT a valid potential function on G4, , p the bound function of 0, P a set of 




U i VPl constraints not implied by cp, and a U i VPl constraint c to be 




added. 




Output: G,^u{c}j its valid potential function tt' and the bound function p' of 




(j) U {c} and the set P' C P of constraints not impUed by U {c}, or 




UiNaAl 11 U |cj- is not satisnable. 


1 


G := G^, TT := TT, p = p, compute -c/(c); 


2 


lor all e t cio 


3 


res . — iNCOONUiFF(^(j , tt , ej, 


4 


11 Tcs — uivo/ii men 


5 


It3LLlIIl 


6 


elsG 


>7 
I 


I tjr , /I 1 . — A ca 


8 


let (it, 1), d) be any edge in E{c); 


9 


compute 5^ and 5^ by using the reduced cost graph for G' via tt'; 


10 


for all X £ V do 


11 


sp := 5^{x) + d + S^{—x); 


12 


p'{x) := min{p(x),LfJ}; 


13 


for all X eV do 


14 


if p'{x) + p'{-x) < then return UNSAT; 


15 


P' := 0; 


16 


for all c' € P do 


17 


{x,y,d') := first element in E{c); 


18 


if 5^ {x) + d + (5^(3/) > d' and 5^{—y) + d + (5^(— x) > d' and 




p'(a;) + p'(-2/) > d' then P' := P' U {c'}; 


19 


return fC?', tt', p', P'; 



Lemma 1 there is equal-weight path from x via the "counter-edge" {—v, —u, d) 
to —X (x~^ = —x~ —V — M ^ — x+ = a;"). 

We calculate the shortest paths in G^ujc} from each vertex x to m 
and from v to each vertex x ((5^(a;)) (line 9). The shortest path for S!^ can be 
computed like 5^ by simply reversing the edges in the graph. 

We can then calculate the shortest path from x to —x via the edge u v 

using the path x+ ^ u w ~^ a;~ as S^{x) + d + x). We update p' if 
required (line 12). 

We can now check satisfiability of </> U {c} in Z using Theorem 3 (lines 13 
and 14). Finally wc check implications using Theorem 4. 

Using the above results, it is not difficult to show that the algorithm is correct 
with the desired complexity bounds. 

Theorem 5. Algorithm 3 (ScSt) is correct and runs in 0{n\ogn + m + p) 
time and 0{n + m + p) space. 



Proof. The algorithm is correct if it returns UNSAT in the case of unsatisfiability 
oi (f> U {c} or the constraint graph G^yjc}; its vahd potential function vr', its 
bounds function p' and the set of constraints P' C P not implied by (/> U {c}. 

The Lemma 2 and the Algorithm IncConDiff (see Cotton and Maler [3]) 
guarantee that after termination of IncConDiff G' = G0u{c} and tt' is its valid 
potential function if ^ U {c} is satisfiable in Q; otherwise (j) U {c} is unsatisfiable 
and the algorithm returns UNSAT. 

After application of IncConDiff the algorithm maintains the bounds func- 
tion (lines 8 to 12) by calculation of the shortest path x u —> v ^ —x via 
one added edge {u,v,d) E E{c) for each node x in G^ujc}- Remark: we only 
have to considered the shortest paths via the added edges p give us the length 
of a shortest path without those added edges. Due to Theorem 3 the algorithm 
checks (f) U {c} for unsatisfiability in Z in the next two lines. If it is unsatisfiable 
ScSt terminates and rctiu-ns UNSAT. 

The remainder of the algorithm computes the set of non-implied constraints 
P' C P by testing for all constraints c' € P if the length of both paths a; u — > 
V y, —y-^u^v-^ —X are longer than d' and the sum of the upper bounds 
p'{x) + p'{—y) is greater than d' where {x,y,d') G E{c'). If all three cases hold 
then c' is not implied by U {c} thanks to Theorem 4. 

The nm-time is determine by the nm-timc of IncConDiff, the calculation 
of (5^, 5^ which are 0{nlogn + m), and the implication check C(p). All the 
other computations can be done in constant or linear time with respect to n and 
m. So the overall run-time is 0{n log n + m+p). □ 

The cost of IncConDiff and the shortest path computations arc each 
0{n log n+m), while the implication checking is 0{p). The space required simply 
stores the graph and implication constraints. 

6 Experimental Results 

We present empirical comparisons of the algorithms discussed herein, first on 
satisfaction and then on implication questions. 

For both experiments we generate 60 UTVPI instances in each problem 
class with the following specifications: the values d range uniformly in from —15 
to 100. approximately 10% arc negative, each variable appears in at least one 
UTVPI constraint, each constraint involves exactly two variables, and there is 
at most one constraint between any two variables two variables are allowed. 

In addition, for the implication benchmarks 10 implication sets P of size p 
were created for each n using the same restrictions as defined above. On average 
over all benchmarks, 65% of the constraint P were implied by the corresponding 

4>. 

The experiments were run on a Sun Fire T2000 running SunOS 5.10 and a 
1 GHz processor. The code was written in C and compiled with gcc 3.2. 

We run incremental satisfiability on a system of m constraints in n variables, 
adding the constraints one at a time. We compare: IncLaMu the incremen- 
talization of LaMu presented in Section 4, ScSt the incremental implication 
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checking algorithm of Section 5 where p = 0, toLaMu running LaMu m times 
for m satisfaction checks, and HaSt the algorithm of [6]. The results are shown 
in Table 2, where d represent the density of a UTVPI instance. We split the ex- 
amples into cases that are feasible, Z infeasible, and Q infeasiblc. Moreover, the 
table entry "all" shows the overall average run-time and the number of exam- 
ples for each case in the same ordering as above written. Interestingly for dense 
satisfiable systems HaSt is best, but overall IncLaMu is the clear winner. 

The incremental implication checked satisfiability and the implications of 
constraints P incrementally as each of the m constraints were added one at 
a time. A run was terminated if there were no more constraints to add, all 
constraints in P were implied, or unsatisfiability was detected. We compare the 
two algorithms that can check implication: ScSt versus HaSt. Table 3 shows the 
results. Overall the checks for implication are cheap compared to the satisfiability 
check for each algorithm. Hence the results are similar to the satisfiability case. 
Again HaSt is superior for dense systems, while ScSt is the clear winner on 
sparse systems. 



Table 3. Average run-time in seconds of the implication algorithms 
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Average run-time of Q infeasible 
problems. 



7 Conclusion 

We have presented new incremental algorithms for UTVPI constraint satisfaction 
and implication checking which improve upon the previous asymptotic complex- 
ity, and perform better in practice for sparse constraint systems. 

Wc can easily adapt the algorithms herein to provide non-incremental im- 
plication checking in 0{n'^ logn -|- nm + p) time and 0{n + m + p) space, and 
generate all implied constraints in 0(n^ logn -|- nm) time and 0{n + m + p) 
space, where p is the number of implied constraints generated. 
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